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ABSTRACT 


Procedures and methods for predicting aerothermodynamic heating to 
delta orbiter shuttle vehicles has been reviewed. A number of approximate 
methods were found to be adequate for large scale parameter studies, but are 
considered inadequate for final design calculations. It is recommended that final 
design calculations be based on a computer code which accounts for non- 
equilibrium chemistry, streamline spreading, entropy swallowing, and turbulence 
It is further recommended that this code be developed with the intent that it 
can be directly coupled with an exact inviscid flow field calculation when the 
latter becomes available. 

The recommended procedure for parameter studies is to calculate local 
pressures based on tangent cone or wedge approximations and heat transfer 
following Eckert's reference enthalpy method. This procedure is relatively 
simple and was found to agree favorably with wind tunnel data for a shuttle 
configuration at angles-of-attack which are of potential interest. 

A nonsimilar, equilibrium chemistry computer code (BLIMP) was used to 
evaluate the effects of entropy swallowing, turbulence, and various three 
dimensional approximations. These solutions were compared with available 
wind tunnel data. It was found from this study that, for wind tunnel conditions, 
the effect of entropy swallo /ing and three dimensionality are small for laminar 
boundary layers but entropy swallowing causes a significant increase in turbulent 
heat transfer. However, it is noted that even small effects (say, ]0-20%) may 
be important for the shuttle reusability concept. 
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SECTION 1 


INTRODUCTION 

The design of an efficient thermal protection system for a shuttle orbi- 
ter requires reliable heat transfer prediction methods. Under orbital entry 
conditions, heating is principally a boundary layer phenomena but even exact 
boundary layer analyses would be inadequate without sufficient methods for pre- 
dicting pressure distributions and edge boundary conditions. These predictions 
of aerodynamic heating and edge conditions ma” be based on analytic methods, 
experimental data or a hybrid combination of these two. 

Analytical techniques, which are presently available, cannot yield pre- 
cise predictions for heating distributions everywhere on the body although, with 
suitable approximations, they can satisfactorily predict heating rates to local 
regions such as the nose, leading edge and centerline portion of the orbiter. 

The key here is "suitable approximations." The test of whether or not an approx 
imate is suitable depends on a comparison of predictions with eith c an exact 
solution or experimental data. An exact solution for shuttle geometries is not 
available and wind tunnel data, though highly desirable, is singularly insuffi- 
cient because of inadequate simultaneous duplication of pertinent parameters in 
ground test facilities and the lack of complete model scaling laws to permit 
confidence in extrapolating to flight conditions. 

Thus the present capability must rely on approximate analytic or semi- 
analytic methods fortified by ground facility experimental data. The ideal kind 
of analysis would be one which is derived from a consideration of all phenomena 
(including for example, homogeneous chemistry) associated with the flight condi- 
tions and which can be generalized to wind tunnel conditions. Then if the analy 
sis adequately predicts the levels of wind tunnel heating, there will be a high 
degree of confidence in the predicted heating to the flight vehicle. 

A very comprehensive review on heat transfer prediction capabilities for 
Apollo-class vehicles is presented in Reference 1. However, shuttle differs 
from Apollo in many aspects, each of which increases the prediction requirements 
The present discussion will center around these difficulties and prediction 
methodology which are directly applicable to shuttle requirements. Two of the 
most outstanding differences between shuttle and Apollo are that shuttle is a 
true three-dimensional body and is large; a typical configuration would cover 



half of a football field. An additional difference which is not readily appar- 
ent is in the design philosophy; each Apollo was flown only once so that it was, 
at least in principle, possible to design early vehicles conservatively and ad- 
just the thermal protection weights on successive vehicles using flight data. 

Such is not the case with shuttle. To be successful, the "first-off" shuttle 
should be designed for upwards of 100 flights so that ultra-conservative design 
procedures would render shuttle impractical. 

The choice of a thermal protection system has a significant bearing on 
the required prediction accuracy. Ablative materials, as were used on Apollo, 
have a high degree of accommodation or tolerance and can survive through sus- 
tained exposure to higher-than-expected heating rates so long as sufficient 
coating thicknesses are provided. In contrast, the surface temperature of a 
radiation cooled surface, as is proposed for shuttle applications, is dependent 
on the heating rates. Metallic surfaces have a maximum service temperatu. e 
which when exceeded results in rapid deterioration of the material and will re- 
duce the degree of vehicle reusability. Hence prediction accuracy is more cri- 
tical for shuttle. 

The selection of materials for the shuttle thermal protection system de- 
pends on the expected peak heating rates which will be maintained lower than 
those experienced by Apollo. This will be achieved by decelerating the vehicle 
at higher altitudes but, at the lower densities associated with these altitudes, 
non-equilibrium chemistry effects are enhanced. These effects may be beneficial 
in the stagnation region if noncatalytic materials are utilized but a penalty 
must be paid in increased downstream heating. Nonequilibrium chemistry also in- 
creases the concentration of dissociated oxygen near the surface which may dras- 
tically increase oxidation rates and hence reduce the life expectancy of the 
vehicle . 

Shuttle size enhances the probability of turbulent flow ii •: ? boundary 
layer which in turn increases heating rates. A serious question aised re- 
garding where on the vehicle and when during the trajectory transition will 
occur. Transition during the peak, heating phase of the ei.try may dictate higher 
temperature materials on the shuttle afterbody. Thus adequate transition cri- 
teria must be specified. 

Apollo was launched in tandam with the launch vehicle and thermally pro- 
tected from ascent heating. Shuttle will probably be launched in a piggy-back 
fashion so that shock interference between the launch vehicle and shuttle may 
locally increase the thermal protection requirements. Further, during entry 
maneuvers, various control surfaces will interact with the normal flow field 
patterns to cause interference heating and flow separation-reattachment problems. 


2 



SECTION 2 


CHARACTER OF FLOW FIELD DURING HEATING 

For purposes of convenience and analytic simplification, assumptions re- 
garding the nature of a hypersonic flow field are often divided into the fluid 
flow and chemistiy aspects. Typical assumptions which are made are shown below. 


Hypersonic Flow Field 


I 

I 

I 




Depending on the entry conditions and location on the vehicle, various combi -.a- 
tions may be simultaneously important. For example, at the stagnation point, 
a laminar boundary layer with nonequilibrium chemistry may occur. The challenge 
of shuttle is that somewhere in it* trajectory there will be separated, inviscxd, 
and viscous flows; the last will involve laminar and turbulent flows and non- 
equilibrium chemical kinetics will be important. In short, all of the above 
phenomena at one time or another will be important. 
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2.1 


STAGNATION REGION 


In the stagnatior region of the nose or leading edge and at low altitudes 
where the density is high, the chemical reaction rates will be rapid enough for 
equilibrium assumptions to be valid. Further it can be shown that there are 
negligibly small errors associated with the assumption that a well defined bound- 
ary layer .xists with zero gradients for edge conditions. Within the confines 
of these assumptions, several solutions (analytic, correlation and experimental) 
are available and have prove accuracy. 

At moderate altitudes, the equilibrium chemistry assumption is not valid 
although the boundary layer assumptions still are. Thus the viscous and invis- 
cid fields are still decoupled and can be analyzed separately. In the inviscid 
shock layer, which is much thicker than the boundary layer, the predominant 
chemical reaction is dissociation. Dissociation, being a two-body reaction, can 
ac.iieve equilibrium even though recombination in the boundary layer is near fro- 
zen. Thus the boundary layer edge conditions can be readily determined without 
solving the inviscid field. Because of this simplification, analyses are avail- 
able for the nonequilibrium stagnation boundary layer. A useful simplification 
and special case of nonequilibrium is a frozen boundary layer. This assumption 
has often been used for studies of surface kinetics. 

At high altitudes, even dissociation rates are not rapid enough to achieve 
equilibrium, so if boundary layer assumptions are retained a more careful analy- 
sis of th ' inviscid flow must be used to. determine the boundary layer edge con- 
ditions. Actually, the boundary layer assumptions commence to breakdown as non- 
equililrium chemistry begins to be important so that at high altitudes nonequi- 
librium viscous shock layer assumptions are used for analysis. Solution methods 
for these flows are available for simplified air models and recently for general 
multicomponent models. 

At still higher altitudes, the shock wave begins to grow and can no longer 
be considered as thin. Eventually the shock wave merges with the shock layer to 
form a viscous merged layer. General chemistry analyses of this flow regime are 
not available but, fortunately, it does not occur during the shuttle heating 
cycle . 


2.2 DOWNSTREAM REGIONS 

Regardless of v‘ 'ther the stagnation flow is of the boundary layer type 
or viscous shock layer type, far enough downstream the windward side of shuttle 
will approach a viscous and an inviscid layer. Because of the chemical reactions 



and effects of shock curvature, the governing differential equations, with 
boundary layer assumptions, are nonsimilar but can still be analyzed provided 
the edge boundary conditions are known. 

However, at high altitudes and just downstream of the stagnation region, 
the boundary layer assumptions may not be valid. For these conditions, the com- 
plete shock layer may be viscous, chemically reacting and nonsimilar, a condi- 
tion for which there are apparently no published exact solutions. 

When the boundary layer assumptions are valid, solutions are available 
for general air, chemically reacting flows around simple body shapes. These 
techniques show promise, with the application of streamline divergence methods, 
of being suitable for shuttle. Nonetheless, even a perfect boundary layer 
analysis would be inadequate if appropriate techniques are not available for 
specifying the edge boundary conditions. Pressure distributions do not seem 
to be very sensitive to chemistry so they are often determined a priori, by 
approximation techniques. Subsequently, the other edge conditions are approxi- 
mated. The ideal approach would be to solve the direct problem for the invis- 
cid field, but present technology* permits this only for simple body shapes and 
even then calculation times are long and not practical for extended trajectory 
studies . 

As the flow progresses downstream, boundary layer growth may lead to 
transition into turbulent flow. Although correlation and analogue methods are 
available for approximate predictions of turbulent heating, the prediction of 
when transition will occur is in a state of flux. The establishment of a tran- 
sition criteria is further complicated by its dependence on the state of the 
gas at the edge of the boundary layer and other boundary layer parameters. 

If the windward side of shuttle poses some difficult prediction prob- 
lems, accurate prediction of the leeward side is virtually impossible. Only 
a very limited number of studies have been conducted on leeward shuttle heat- 
ing. At very small angles-of-attack , when the flow remains attached, methods 
similar to those Used for the windward side are applicable, although it should 
be noted that here the streamlines are converging rather than diverging as on 
the windward side. However, at high angles of attack, flow separation occurs and 
the viscous-inviscid interactions are not yet amenable to analysis. Extensive 
experimental data and conservative extrapolations appear to be the only recourse. 

•ff 

In this context, "present technology" means existing, operational and reliable 
computer codes. 
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SECTION 3 


APPROXIMATE HEAT TRANSFER PREDICTION METHODS 

3.1 STAGNATION AND LEADING EDGE REGIONS 

The state-of-the-art for predicting stagnation region heat transfer 

rates to planar or axisymmetric bodies is well m hand. Adequate predictions 

2-11 

can be made using available solutions for both non-catalytic and fully 

catalytic surfaces so that only a cursory discussion will be presented. 

Assuming a binary air model, reference 2 also presents general nonequilibrium 
solutions for surfaces of arbitrary catalycity. Two difficulties occur when 
arbitrary catalycities are concerned. First, there is very little data on 
the dependence of surface catalycity on temperature and second, there is no 
information to adequately define catalytic efficiency for gases with more 
than one dissociated specie. However, an important conclusion to be reached 
from examination of these analyses is that for stagnation point flows to sur- 
faces of infinite catalycity at temperatures of interest there is only a 
negligible difference in heat transfer between assumptions of equilibrium and 
nonequilibrium chemistry. This conclusion is especially important because of 
the simple correlation solutions which have been obtained for equilibrium flows. 
Thus, since experimental evidence shows chat metallic surfaces tend to be cata- 
lytic, these equilibrium chemistry solutions can be used to predict the heat 
transfer rates. In addition, for low catalycity surfaces, these predictions 
represent an upper bound on the expected heating rates. 

The accuracy of the above analyses have been adequately proven for both 

ground test and flight hardware, but in order to apply them to shuttle whose 

stagnation point may be neither axisymmetric nor planer, some method such as 

reference 12 must be used to account for three dimensional effects. For in- 

2 

stance, for an axxsymmetrxc stagnatxon poxnt xn a dissociating flo* 



where 


n 


0.52 equilibrium chemistry 
0.b3 frozen chemistry 
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Then for R.^, and R c as the principal radii of curvature of a non-axisymmetric 
stagnation point. Equation (3-1) is scaled by the expression 
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and [Nu x //Re x ] 

equal to A Note that for R 


T* 

stagnation line. 


is the stagnation point value for 
>> R t the solution 


an axisymmetric nose radius 
degenerates to that for a 


The leading edges of the delta wing portion of the vehicle are not stag- 
nation lines and should be considered as swept or yawed cylinders. Solutions 
for finite length cylinders are not available since this comprises a true three- 
dimensional bod/. In addition, nonequilibrium chemistry has not been generally 
considered. Even so, for catalytic surfaces, adequate predictions can be made 
by assuming that the leading edge is an infinite swept cylinder. Then, draw- 
ing upon stagnation point experience an equilibrium chemistry correlation solu- 
tion such as Reference 6 can be used. This will be discussed in Section 3.2.2. 


Overall, it is generally concluded that, with the exception of surfaces, 
of finite catalycities , the confidence level is quite high, even for viscous 
shock layer flows, for the prediction of stagnation and leading edge heat 
transfer. 


3.2 DOWNSTREAM REGIONS 

Although it may be necessary to consider viscous shock layers in the 
vicinity of the stagnation point, the major portion of the lower surface of 
shuttle can be characterizes as having an inviscid region and a boundary 
layer region. Under these conditions it is the boundary layer that is of 
primary interest, since an adequate analysis of it would provide required 
design heating rates. However, the inviscid flow, either in a coupled or 
independent sense, is important because it provides the boundary conditions 
for the boundary layer equations . 
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3.2.1 INVISCID FLOW 

There are presently some solutions 13 " 20 for the inviscid flow field 

about symmetric bodies at angles of attack but a general solution for an arbi- 

21 22 

trary three-dimensional body is not readily available, although some methods ' 
are in the development state. For a parametric study of entry trajectories, 
the inviscid flow field or boundary layer edge conditions should be determined 
by approximate means; more exact methods being reserved for final trajectory 
calculations. With the inviscid flow field and the boundary layer assumed to 
be uncoupled, approximate techniques can be used to specify the boundary layer 
edge conditions and a solution of the inviscid field becomes unnecessary. 

23 24 

It has been noted ' that pressure distributions are relatively 
insensitive to chemistry and displacement effects at sufficiently high 
Reynolds numbers so that most approximation methods center first upon obtain- 
ing the surface pressure distribution and then using it to determine the chemi- 
cal, thermodynamic and dynamic states of the fluid at the edge of the boundary 

25 

layer. For simple axisymmetric bodies, it has been shown that a modified 

Newtonian pressure distribution is sufficient (within 5% of experimental data) 
and is the assumption used by several investigators 10,11 * 25,26 . However 

27 

for wind tunnel delta winged models, the modified Newtonian pressure was found 
to be 15 to 20% lower than experimental measurements. For this configura- 
tion, Marvin et al show that the elliptic-cone technique is more accurate. 
However, Fannelop 30 shows that a similar but simpler procedure (effective cone 
technique) is adequate for symmetric bodies at angle of attack by comparing 
predictions with the experimental data of Cleary 31 . Laminar heat transfer 
rates are approximatley proportional to the square root of the pressure so that 
relative errors in heat transfer rates, due to inaccurate pressures, will be 
less than the relative pressure errors. Thus, pressure approximations which 
are accurate to the order of a few percent should be sufficient and more 
refined approximations would be of second order and would probably be over- 
shadowed by the question of boundary layer edge chemistry. An exception to 
the above would be for flight conditions for which turbulent flow is expected. 
Transition phenomena is dependent on the edge Mach number which in turn depends 
on the pressure distribution. Experimentally determined boundary layer edge 
Mach numbers on a straight-winged orbiter model were compared with those pre- 
dicted from a tangent cone pressure distribution in reference 32 and found to 
be in good agreement. It then appears that a tangent or effective cone approx- 
imation for pressures should be sufficient for shuttle centerline predictions. 
There is insufficient data to reach any conclusions about pressure distribu- 
tions on the winged portion of the orbiter although here again tangent cone 
or, more probably, tangent wedge approximations should be adequate. 



Given a pressure distribution on the vehicle centerline, one possible 
approach for determining the edge conditions would be to calculate an isentropic 
expansion of the flow from the stagnation point to the local body pressure. The 
chemistry may be considered as being in equilibrium, nonequilibrium or frozen, 
with the choice depending on flight conditions. For a significant portion of 
the shuttle trajectory the altitudes will be high enough so that recombination 
rates will be relatively slow, hence the expansion will be in chemical nonequi- 
librium. Even so, an equilibrium expansion is often used since, for catalytic 
surfaces, it is known that stagnation point heating is relatively insensitive 
to chemistry, and the same is presumed to be true elsewhere on the body. Fur- 
thermore, equilibrium calculations are simpler and edge chemistry does not di- 

5 33-39 

rectly enter into some of the heat transfer prediction schemes. ' But, 

for r.oncatalytic surfaces, nonequilibrium expansions for sphere cones at zero 

40 

incidence were suggested by Blottner and were shown to be of significant im- 

25 

portance by Lewis for a hyperboloid under flight conditions which simulate 
points on a typical shuttle trajectory. 

Even a nonequilibrium adiabatic expansion may be inadequate for shuttle 

because of its largo size. Far downstream of the stagnation point, the bound- 

dary layer edge gas does not originate from the normal shock; rather, it 

originates from an oblique shock. A similar situation arises for flows with 

large shock curvature since an expansion, of any kind, from the stagnation 

region presumes that the entropy at the edge of the boundary layer has the 

high value generated by a normal shock wave. In reality, at a sufficiently 

far downstream location the stagnation entropy will be gradually swallowed 

by the boundary layer and the edge entropy will decrease and approach the 

41 42 

value behind an oblique shock . Hamilton compares the assumptions of an 
edge condition determined by expansion from the stagnation point and the state 
determined by assuming that the flow has just passed through an oblique shock 
of sufficient strength to get the local pressure. As might be expected, the 
lower edge entropy for the oblique shock assumption was shown to cause an 
increase in the predicted heating rates. 

Because of the effects induced by shock curvature and entropy swallow- 
ing, a complete vorticity interaction analysis would be highly desirable, 

albeit very difficult. Because of these difficulties, some semi-coupled 

43 

approximate techniques have been developed. For example, Adams determined 

the shock shape and pressure distribution in the inviscid flow using an equi- 

44 

librium chemistry model and the methods described by Lomax and Inouye . Then, 
assuming that pressures and shock shape are insensitive to chemistry, the non- 
equilibrium state was determined by integration along streamtubes using this 
predetermined pressure distribution. Finally, boundary layer solutions were 
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obtained by iteration of the mass flow within the boundary layer and the 

45 

absorption of inviscid streamtubes as per Kaplan . It becomes apparent 
that, even approximate coupling techniques are very complicated and probably 
very time consuming and application to more complex bodies is hampered still 
further by the need for an adequate description of a three-dimensional entropy 
layer. 

A more exact accounting of shock curvature and nonequilibrium chemistry 

effects requires very extensive computer codes which are not available for 

general three-dimensional bodies. However, codes are available for bodies in 

14 17 

chemically equilibrium flows ' but even under these conditions CDC 6600 
computational times for sphere-cones at angle of attack are in excess of 
5 minutes.^® For nonequilibrium chemistry and small angles of attack., a per- 
turbation technique which uses available zero incidence solution techniques 
such as that used by Fannelop^ could be employed. In general, unless ade- 
quate similarity laws are developed, any technique which uses "exact" solutions 
of the inviscid flow would require a large computer expenditure. 

For small angles of attack and far downstream of the stagnation region, 
the boundary layer edge condition is largely determined by the state behind a 
weak, oblique shock. Then the degree of dissociation in the inviscid flow will 
be law ar.d ideal gas assumptions , along with a specified pressure, may be suf- 
ficient for determining the edge state. This, of course, does not preclude 
large amounts of dissociation caused by viscous heating within the boundary 
layer. 


3.2.2 BOUNDARY LAYER 

With a specified boundary layer edge condition, several approximate 
methods may be used to predict the heat transfer rates. These include Eckert's 
reference enthalpy ^ ' ; Reynolds analogy^ rho-mu^' 50 ; or swept and modi- 

fied swept cylinder methods and are representative of methods developed by a 
combination of theory, experiments and intuition. 

Eckert's reference enthalpy method is based on the assumption that 
incompressible constant property solutions can be used to calculate compress- 
ible flow heat transfer if properties are evaluated at an appropriate reference 
temperature or reference enthalpy. This reference enthalpy is defined as 

i* - i e + 0.5<i w - i e ) + 0.22 (i r - i e ) 

At a stagnation point the recovery enthalpy i r is approximately equal to the 
edge enthalpy i so that 
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For other regions 


i* « 0.5(i e + i„) 


u e 

*r “ L e * 1 T 


For laminar flows r ^ /pr and for turbulent flows r * y'Pr. Typical imcompressible 
constant property solutions are j 
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(stagnation point) 
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(stagnation line) 


xn 



0.332 (Pr) l/? 



0.185 Re Pr 

A 


l/3 


(1 °9l0 Re x> 


2.584 


(laminar flat plate) 


(turbulent flat plate) 
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Re = 
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p e u e x 


In the reference enthalpy method these are then written as 


Nu* 

— — = 0.763 (Pr*) 0 * 4 (stagnation point) 

/rSJ 

NU* 

— — = 0.570 (Pr*) 0 * ** (stagnation line) 

/SSJ 


Nu* . 

— — - = 0.332 (Pr *) l ' 3 (laminar flat plate) 

/KeJ 


0.185 Re*Pr*‘/ ! 

Nu* = (turbulent flat plate) 

x (log^Re^) 2 - 58 * 
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For dissociated air, the Prandtl number is nearly constant so that the reference 
enthalpy equations can Le written in an alternate form by noting that 



1/2 


Using the experience of the analytic studies performed in Reference 2, the lami- 
nar solutions can be corrected for dissociated gases and non-unity Lewis numbers 
by multiplying the right hand sides by the function 


1 


(Lr 11 



where 


n 


0.52 equilibrium chemistry 
0.63 frozen chemistry 
| 0.50 couette flow 


Then, by approximating any point on the body as a stagnation region or a zero 
pressure gradient wedge/cone, the local heat transfer can be calculated. 

The Reynolds analogy method assumes that the velocity and temperature 
fields are proportional to obtain a simple relationship between heat transfer 
and skin friction. The Reynolds analogy factor is then defined as 


K = 


2 

C 


St 

f 


The successful use of the Reynolds analogy depends on an adequate solution for 
the skin friction coefficient and a method of calculating the Reynolds analogy 
factor K. Some typical values of K are 


K = 1.0 


Assumes Pr = 1, M z 0 


*K = 



c;’ 


Empirical determination by Colburn"' 1 


1 + 


K = 


r + 


vv 

Pr(u es /u e> 


o 



Theoretical value determined by 
Rubesin^ 8 
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exp 



1.561 (H - ^M’) 


Empirical determination by 
Reshotko and Tucker^? 


- 1.3 



K 



Von Kdrmcin 


For laminar flow over a flat plate with constant properties, the skin friction 
coefficient is 


C 


f 


0.664 



For turbulent flows C f may be determined with methods set forth by Spalding and 
35 1 53 34 

Chi, Sommer and Short, Van Driest or may use methods such as rho-mu or 

Eckert's reference enthalpy. The number of possible predictions obtained from 

a permutation of K and C f is very large. Select combinations have been com- 

1 54 

pared with experimental data by Pearce for flat plates and by Hopkins and 
55 

Inouye for flat plates and cones. Pearce concludes that the Spalding-Chi 
correlation with von Kcirm^n's Reynolds analogy factor is best, whereas Hopkins 
and Inouye conclude that the Van Driest analysis with a Reynolds analogy factor 
of 1.0 is best. These differing conclusions are typical for turbulent flows 
and are indicative of the fact that turbulence is not a well understood phe- 
nomena and the results are highly dependent on the particular application. 

The swept cylinder methods are based on the analyses of Beckwith,** Beck- 

OO Cf 

with and Cohen, and Beckwith and Gallagher. An important conclusion drawn 
from these solutions is that, even for relatively large sweep angles, the ratio 
of local heat transfer coefficient to leading edge heat transfer coefficient 
along a plane perpendicular to the leading edge of a yawed circular cylinder 
is insensitive to the angle of yaw. For wall temperatures at which there is 
no dissociation the ratio of heat transfer coefficients based on enthalpy 
be expressed as 



can 



where 


t 


i + 


u 2 


+ w 2 
2 


2 *s du e 
our 3 


/ p w> J w u e r2:ldlt 


and 0 * and 0 ' can be determined from tabulated results presented in Reference 
W W f s 

33. The local value of heat transfer is then calculated from 


<*w <*w (i Aw *w^ 


with 


Vs - 0-577 <Pr w >-"-‘ (p w v w ) 


0 s 


‘Vs’ 


0.44 



5 


and 


u_ + v' 


Aw 


= i + r 


The rho-mu method^® is based on integral solutions of the momentum and 
energy equations and uses boundary layer thickness parameters and a reference 
density-viscosity as variable functions. These functions were determined from 
available exact similarity solutions and, with suitable modifications, account 
for cross-flow gradients, streamwise gradients, nose bluntness and real gas 
effects. Hand calculation of the sot of rho-mu equations, though possible, 
is not practical. To facilitate hand calculations, simplified correlation 
appro. ..nations are presented in References 50 and 57 for plates, cones, swept 
cylinders and the centerlines of sharp delta. In addition, modification pro- 
cedures for streamline divergence and variable wall temperature cases were 
described. 
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SECTION 4 

CURRENT METHODOLOGY 

Current state-of-the-art for predicting shuttle heating rates is not up 
to the task of considering the vehicle as an entity because of it? complex 
shape. In lieu of this capability, the vehicle is usually segmented into sec- 
tions for which the heating rates are separately analyzed or approximated. 

Typical surface regions that are considered are shown in Figure 1. The stagna- 
tion region, leading edges and windward surfaces have received the most atten- 
tion since these regions are expected to be subjected to the highest heating 
rates and, fortunately, are easier to approximate than the leeward side. How- 
ever, due to an interaction between separation-reattachment phenomena and tran- 
sition to turbulence, certain portions of the lee side have been noted to have 
significant heating rates. 

4 . 1 WINDWARD SURFACE 

Even though Marvin et. al. show that Newtonian approximations are inade- 
quate for the NAR delta-wing model, the approximation is often favored because 
it can be expressed in a convenient an*.. . orm. Moreover, Newtonian approxi- 

mations are known to be adequate for sphere-cones and other simple blunt bodies 
and are apparently also adequate for lift and drag calculations of delta vehi- 
cles. However, one would expect that, since Newtonian approximations (based 
on body angle) predict the same pressures for cones and wedges, they would be 

inadequate for regions far from the stagnation point. Hence, depending on the 

59 

region of interest, different assumptions may be required. Young et. al. , 
believing that transition to turbulent flow is enhanced by entropy layer swallow- 
ing, use a high entropy Newtonian approximation in the laminar flow region and 
the lower entropy tangent wedge approximation in the turbulent flow region. A 
tangent wedge approximation was also used in Ref ere ice 60 but boundary layer 
edge approximations were modified with an empirical gas constant. Guard and 

Schultz 6 * based their pressures on a combination of experimental data and blast 

62 63 

wave theory of Creager with a correction for nose bluntness. Thomas takes a 

novel approach; he uses a heat transfer prediction technique and model center- 
line heating data and works backwards to obtain the pressure distribution which 
yields the best correlation. With this technique he shows that blast wave cor- 
relations with nose bluntness effects are best at slightly negative angles of 
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attack, irodified oblique shocks are best at zero incidence and a modified 
oblique shock accounting for streamline divergence is best at slight positive 
angles of attack. 

Although approximate pi diction methods have shown generally good agree- 
ment with low enthalpy wind tunnel data, questions dealing with the effects of 
chemistry, shock curvature and flow interactions need yet to be adequately re- 
solved for flight vehicle prediction confidence. Even for wind tunnel models, 
pressures along the windward centerline are the most often measured so that more 
effort needs to be expended on pressure measurements on the remainder of the 
vehicle. 

Given a pressure distribution, an isentropic expansion or oblique shock 
conditions are the simplest methods for approximating the boundary layer edge 
flow conditions. However, in efforts to obtain improved approximations, some 
investigators have employed more sophisticated methods to account for three- 
dimensional effects. 

Marvin et.al. assume that the edge conditions on the vehicle centerline 
are the same as those behind a swept cylinder inclined at an angle equal to the 
local body incidence. However, on the wing portion of the vehicle the edge con- 
ditions are assumed to be determined by an isentropic expansion from the leading 
edge. Young et.al. use an isentropic expansion for low angles of attack whereas 

at high angles of attack they use a swept cylinder theory but correct the stag- 

64 

nation line velocity gradient to account for noncircular cross sections. Masek 
54 

and Pearce account for the stagnation line cross flow by assuming that the 
cross-flow velocity gradient is equal to that which occurs on a circular disk 
with a radius equal to the local wing semi-span and a normal velocity equal to 
the local normal component of the free stream flow. Pearce points out that 
the cross flow correction to a plane oblique shock varies from 0 to 13 percent 
for incidences between 0 and 50° which, at least for high angles of attack, 
makes the correction significant. The data and heating predictions of Marvin 
et.al and Guard and Schultz represent most of the reported windward side, off 
centerline heating data for delta vehicles so that not much can be reported on 
methods applied to the wing portion for determining boundary layer edge condi- 
tions. Whereas, Marvin et.al. predict a three-dimensional edge condition, as 
noted above. Guard and Schultz do not. Guard and Schultz use a planar oblique 
shock for edge conditions and correct the heat transfer results to account for 
cross flow. 

Current methods for predicting windward-centerline, laminar heat transfer 
rates arc mostly based on adaptations of correlation methods such as Eckert's 
reference enthalpy or rho-mu theories. Moote,^ for instance, used the 
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reference enthalpy method but modified it for conical flows. At low angles of 
attack Young et.al. also used the reference enthalpy methods but at high angles 
of attack, a swept cylinder theory (which is a spherical stagnation point theory 
modified for two dimensionality and sweep) was used. Marvin et.al. used the 
Beckwith and Cohen cross flow theory to modify the finite difference scheme of 
Marvin and Sheaffer 6 * 5 for streamline divergence effects. This was found to be 
in good agreement with data for regions aft of the wing-body junction but for- 
ward of this point, swept cylinder theory was found to represent the data better. 

For laminar flow on the wing of the vehicle, the cross flow theory of 
Beckwith, and Beckwith and Cohen were used by Marvin et.al. and Young et.al. 

For lifting body configurations. Guard and Schultr used a similar approach for 
the underside heating whereas Reference 60 uses Eckert's reference enthalpy 
method and flat plate solutions with a S’.rip theory correction for cross flow. 

42 

For turbulent flow along the vehicle centerline, Hamilton used the 
Eckert reference enthalpy method with an origin for turbulent flow which is 
assumed tc occur at the s Lai L of Lran^i Liuu . The Spalding -Chi method was used 
for simulated shuttle vehicles at low angles of attack by Young et al., Masek 
and Forney and Moote. The last modified the method used to account for coni- 
cal flow and real gas effects by using an empirically determined weighting fac- 
tor of 1.25. At high angles of attack Young et al. chose to use the results of 
Beckwith and Gallagher, whereas Marvin et al. elected not to be selective and 
found good agreement between their data and the theories of Spalding-Chi , Sommer- 
Short and Van Driest for Reynolds analogy factors of 1.0. Hopkins and Inouye 
recommend Van Driest' s method with a factor of 1.0 but note that the Spalding- 
Chi method with a factor of 1.2 also correlates well with data. 

Turbulent flow on the leading edge was predicted by Guard and Schultz 

by representing it as an isolated swept cylinder and applying the results of 

Beckwith and Gallagher. Turbulent heating throughout the windward side of a 

lifting body was predicted by Reference 60 using Eckert's reference enthalpy 

68 

method with a modified Reynolds analogy and the Schultz -Grunow skin friction 

law. 

It is expected in both laminar and turbulent flows, that the accuracy of 
the predicted heating rates would be dependent on the accuracy of the specifica- 
tion of the boundary layer edge condition but Young et.al. and Pearce point out 
that the transition criteria is also dependent on the edge condition. Thus wind 
tunnel free stream conditions may cause transition to occur earlier on models 
than might be expected in free flight so that conservative estimates would be 
obtained by the direct use of the tunnel data. 
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4.2 

For very small angles of attack, where the flow on the leeward surface 
does not separate, the heating predictions have used procedures similar to chose 
used on the windward surface. For example. Guard and Schultz predict the leeward 
heating for a delta-body vehicle using two-dimensional rho-mu theory and noted 
good agreement for both laminar and turbulent flows at angles of attack up to 
30°. The good correlations suggest that separation-reattachment phenomena* does 
not occur in the vicinity where measurements were made. However, the data of 
Maise 6 ^ for semi-pyramidal shapes at angles-of-attack less than 55° and the 
data of Hefner and Whitehead for a delta-wing orbiter at angles of 20° and 40° 
show that the leeward side heating may be nonuniform with regions of "peak" 
heating. This nonuniform heating is apparently a result of three-dimensional 
separation-reattachment phenomena complicated by transitional and turbulent 
flows. Under these conditions it is generally conceded that the prediction of 
this leeward heating, which requires estimates of pressure and edge flow condi- 
tions,** is much too ambitious for current and even near-future technologies. 

As a first approximation. Reference 60 suggests using tangent cone approxima- 
tions when the local inclination is positive and use flat plate theory for re- 
gions of zero or negative incidence. In this way, transition and turbulent ef- 
fects can be included but separation-reattachment possibility would cast some 
doubt on the validity of the predictions. 

The only recourse to the above dilemma appears to be semi-empirical 
techniques utilizing wind tunnel data. Although this would cause a decrease 
in design confidence, the leeward heating is generally sufficiently low so as 
not to be serious. Even so, it would not be desirable to be overconservative 
in the selection of surface material since the leeward area is at least 50 per- 
cent of the total vehicle area. 

4.3 SHORTCOMINGS OF CURRENT METHODOLOGY 

With the exception of the stagnation point, the current prediction meth- 
ods are founded on modifications of solutions which do not, in a physical or 
mathematical sense, completely account for the real environment to which shuttle 
will be exposed. This, of course, does not imply that these procedures have 
been ihadequate for preliminary design purposes since a large number of conditions 
__ 

Which is one of the reasons why the authors chose the delta-body rather than 
a delta-wing vehicle. 

** 

In passing it should be noted that even the prediction of flow fields for 
simple shapes such as sphere-cones at angle of attack is no simple task and 
is the objective of much activity. 69-72 
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were studied. However, for final design purposes, the present methods have 
certain shortcomings that must be minimized or the methods must be replaced by 
more adequate analyses. A complete inviscid/viscous flow field analysis would be 
ideal but such a program is not available and its development may not be possible 
with the resources allotted to shuttle development. 

As already discussed, the prediction procedure for the windward side can 
be considered in three different parts. First, the pressure distribution is de- 
termined, then the boundary layer edge conditions are calculated, and finally an 
appropriate boundary layer solution is used to predict the heating rates. Methods 
for the prediction of surface pressures are believed, on the whole, to be suffi- 
ciently accurate. The other two parts are inadequate for certain critical por- 
tions of the shuttle trajectory. In establishing the edge conditions, the two 
most common assumptions are that the flow originated either from a stagnation 

region or an oblique shock. The former is valid near and the latter far from 

43 

the stagnation region. Adams shows that "far" may be as high as 15 nose radii 
dovnstream for sphere-cone bodies. In between these extremes, say between 3 and 
15 nose radii, a three-dimensional entropy layer will form and swallowing by the 
boundary layer may be important. None of the above methods accounts for this 
entropy swallowing effect. In addition, present methods rely on solutions ob- 
tained with similar and locally similar assumptions when in fact the boundary 
layer will be highly nonsimilar. 

A somewhat more serious question arises from flow field chemistry. The 
present methodology makes use of analyses developed for ideal gas flows and cor- 
rects for real gas effects by empirical methods (e.g., Eckert's reference en- 
thalpy method) . Although these precedures have been shown to yield good agree- 
ment for stagnation flows with equilibrium chemistry, their application to down- 
stream flows of complex bodies is not fully justified and in nonequilibrium 
chemistry with noncatalytic surfaces, these methods fail except to bracket the 
probable heating range. Moreover the prediction methods are based on similarity 
or local similarity assumptions which are generally invalid for nonequilibrium 
chemistry boundary layers. 

Turbulent boundary layer heat transfer predictions may be very important 
for shuttle. As noted though, there are several prospectively good approxima- 
tion methods; the choice of which may depend on the particular geometry and/or 
environmental conditions. None of these methods is believed to be universally 
correct; in fact, since none of these methods directly account for nonsimila~ity 
and thermochemistry effect , they i«'.ay be seriously in error for shuttle even 
though they adequately predict wind tunnel heating rates. 
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RECOMMENDED CURRENT STATE-OF-THE-ART METHODOLOGY 
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It has been noted that pressure distributions are relatively insensitive 
to chemistry so that pressures can be determined from ideal gas flows in the 
form of experimental correlations or analyses. Modified Newtonian flow methods 
are the simplest to use since the pressure can be expressed analytically, but 
it was shown in wind tunnel measurements on a scale delta-wing vehicle that the 
predicted centerline pressures are about 15-20 percent lower than measured val- 
ues. Nonetheless, for rough-cut calculations on the vehicle centerline, this 
should be sufficient. For more accurate predictions a tangent cone method is 
recommended over the more complex elliptic cone method. A comparison of these 
two predictions and the data of Marvin et.al. is shown in Figure 2 and good 
agreement is obtained up to an angle-of-attack of 40°. At 53.5° the tangent 
cone method overpredicts the pressure by about 10 percent whereas the elliptic 
cone method underpredicts by about 5 percent. Thus the tangent cone approxima- 
tion will lead to conservative (higher) heat transfer predictions. In retro- 
spect, the failure of the tangent cone approximation at 53.5° should be expected 
since this angle approaches the maximum cone angle that will support an attached 
shock and, at these high angles, the precedure would begin to fail even for 
slightly blunted cones. 

If an approximate analytic equation is desired, a modified tangent cone/ 
wedge method can be developed as follows. The Newtonian pressure is given as 

p - p = p U 2 sin 2 ct 

where a is the angle of the surface with respect to the free stream vector. At 
a stagnation point the Newtonian pressure is higher than the actual pressure 
whereas on sharp wedges and sharp cones it is lower. The modified Newtonian 
pressure adjusts the predicted value to be correct at the stagnation point; 
then the pressure is given by 


P - Pc. 
P s “ Po 


= sin a 


This modified Newtonian pressures does a very good job of predicting the pres- 
sure distribution on the forward portion of blunt bodies, but far from the nose 
region this procedure predicts a slightly lower pressure than the unmodified 
method. Thus the discrepancy is widened further. Since the Newtonian formula- 
tion is analytically convenient the following procedure is recommended for large 
ratios of surface distance to nose radius. Let 
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P ' Poo 
P r - Poo 


sin 2 a 

sin 2 a r 


where p r is the pressure on a reference surface oriented at the angle a^. This 
equation can be written in the form 


E_ = 


P 


00 



sin 2 a 

sin 2 g r 


( 4 - 1 ) 


which is then a Newtonian tangent cone/wedge equation. To illustrate the use 
of this equation, consider an airfoil at angle of attack as shown below 



The chord line will be used as a r and, depending on the nature of the flow around 
the airfoil and whatever body is attached to it, the pressure p r can be determined 
for a wedge or cone at angle g r> The pressure at A which has a local angle of 
incidence equal to a can then be calculated from Equation (4-1) . For a typical 
shuttle vehicle this procedure will yield pressures comparable to tangent wedge 
or tangent cone values as shown in Figure 3. Note, however, that a qualitative 
decision must still be made regarding which, tangent cone or tangent wedge, 
is most applicable. Note also that the procedure would not be valid in the 
limit as a r approaches zero. 

At low angles-of-attack the transverse curvature of the shock wave on 
the windward side will not be significant except near the centerline of the vehicle. 
The flow across the wing is then akin to that around a yawed-blunted wedge so 
that the pressures should be adequately predicted with tangent wedge approximations. 
At higher angles-of-attack the shock wave, even over the wing, will have signifi- 
cant transverse curvature so that tangent cone approximations should be used. 

It is apparent then, that a decision must be made as to when each assumption is 
valid. Based on the Marvin data, the following rule-of-thumb is recommended: 

If the angle-of -attack is somewhat greater than the half-angle of the delta, 
then use the tangent cone approximation; and if the angle-of-attack is approxi- 
mately equal to or less than the half-angle, use the tangent wedge approximation. 
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Regardless of which approximation is used on the wing, the cor oination 
of the airfoil shape and dihedral causes the true local angle of incidence to 
be different from the vehicle angle-of -attack. For small values of the sum 
(a + a) , a geometric analysis will yield the following results for chordwise 
stations greater than a few percent. 

sin y = cos <p sin (a + a + a) 


where 

a = vehicle reference line of attack 

a = angle between chord line and reference line 

a = local surface inclination with respect to chord line 

<j> = local wing dihedral angle 

y = true angle of incidence on wing 


Calculations using y in the tangent wedge approximations are compared with the 
Marvin data in Figure 4. It should be noted that the general decrease in prep- 
sure going outboard is due primarily to a change in a rather than a three-dimen- 
sional end flow effect. 

The fins can be handled in much the same way as the wings where once 
again the true incidence of the fin must be determined. Assuming that the fin 
is a flat plate with no airfoil shape, this incidence is given by 


sin y 1 “ cos sin 2 a + tan 2 <j>* sin(£ + a) 


where 


a 

Y' 


sin £ 


sin a 

^/sin 2 a + tan 2 4>'' 


angle between *in chord line and vehicle symmetry plane 
angle of tilt of fin 
true angle of incidence 


There is no available data on fin pressure distributions to use for comparisons 
of predicted values. 

The pressure along the stagnation line of both the fin and the wings 
are calculated as the pressure on a swept cylinder with a correction for angle- 







of-attack to obtain the true yaw angie. That is, for the wing, with a small 
dihedral angle 


where 


cos A = cos a cos C 


5 *= semi-apex of wing 

A « true angle of wing leading ledge 

and for the leading edge of the fins 

cos A' = r sin (6 + a) 


where 

r = ( 7 + (sin <5 tan a + cos 6 tan <J> ' 

6 = nominal sweep angle of fin measured on vehicle symmetry plane 
A 1 = true angle of attack of fin leading edge 

The above angles fcr the wing and fin are shown in Figure 5 and the above rela- 
tionships are derived in Appendix 1. 

As far as the leeward side is concerned, short of a complete flow field 

analysis, there are no adequate means of predicting surface pressures. Based 
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on the analysis of Reference 73 and experimental data, it is recommended 
that laminar heating rates be calculated as 0.56 of the laminar flat plate 
value and turbulent rates at 0.84 of the turbulent flat plate value. 

For the lovr velocity-low altitude portion of the shuttle trajectory, 
the predicted pressure can be used with equilibrium gas assumptions to obtain 
adequate boundary layer edge conditions. But during the earlier portions of 
the trajectory, accurate predictions of the edge condition must include the 
effects of nonequilibrium chemistry, especially if surface kinetics are to 
be considered. However, for a "f x *st cut" design calculation, equilibrium 
chemistry can be assumed with the knowledge that the use of noncatalytic surfaces 
affect the heating rates both at the non-catalytic surface and downstream of it. 
For a distance up to several nose radii downstream, the edge condition can be 
determined by a streamline expansion. 
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du 
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e dx 


1 dp e 

3bt 


it = ie + ~2~ 


= constant 


ie ie ^ p e ,? e^ 
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Far downstream (say a distance in the order of tens of nose radii) the edge 

condition can be assumed to be the same as that which would exist behind an 

oblique shock with p 2 equal to the local surface pressure. Thus with p 2 /p 1 

known, equilibrium normal shock tables can be used to determine the equivalent 

free stream Mach number, M , which is equal to M sin 8. The tables will 

equiv. 00 

also yield all static properties and the normal component of velocity; thus 
with the tangential terms unchanged the edge velocity can be calculated to 
complete the edge conditions. For intermediate distances, there will be a 
transition from the high entropy stagnation state to the low entropy oblique 
shock state. The rate of transition depends on the shock wave shape and is 
further complicated by three dimensional effects. Since the heat transfer 
prediction methods to be recommended use similarity or local similarity assump- 
tions, the results are not sensitive to the smoothness of the edge conditions. 
Thus it is suggested that, for moderate distances from the stagnation point, 
the solution be bracketed by solving both the normal shock expansion and the 
oblique shock conditions. 

It is believed that the particular choice of an approximation method 
to be used to predict the heat transfer rates is not critical since a 
degree of empiricism is built into each method. In line with the tangent 
cone/wedge method used to predict pressures, it is recommended that the heat 
transfer be predicted using, for example, Eckert's reference enthalpy method 
for conical or wedge flow where the local body incidence is used as the cone 
or wedge angle. This procedure, of course, would not be valid near the nose 
or leading edge but would be applicable at large S/R. Calculations, as described 
are compared with the Marvin's centerline data for a = 15° and a = 30° in 
Figures 6 and 7 respectively. The 53.5° case is not compared since, as noted, 
this is too close to the limiting cone angle. The agreement at a - 30° is very 
good for X/L 0.2 but at a = 15° the predicted values are conservatively higher 
than measured values. 

Stagnation point heat transfer can also be calculated using the reference 

enthalpy method but between the stagnation point and x/L ^ 0.2 the method of 
76 

Lees is suggested. Briefly, the ratio of local heat tiansfer to stagnation 
heat transfer is given as 

q _ F(s) 
q o ~ F (oY 


I . 

e.g,, Reference 75. 
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two-dimensional 
axi symmetric 


and the coordinates are defined by the following sketch: 
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An alternative approach for regions not too close to the nose is that 

used by Guard and Schultz. On the windward centerline this procedure calculates 

the two-dimensional heating rate and scales it with methods developed by 
77 57 

Dunavant and Thomas, et. al. For example, the local heating can be predicted 
by Eckert's method for flow over a plate or wedge, i.e.. 


— - = 0.332 (Pr*) 
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then the ratio of heat trasnfer coefficients is 


where 



2 (tan eXM^cos a) 


(4-2) 


(4-3) 


and e is the half angle of the delta wing. The above procedure would be valid 
only for the centerline of the delta portion of the vehicle. When turbulent 
flow is expected, it is observed from Equation (4-3) that the effect of stream- 
line divergence is very small and can probably be neglected. 

For some vehicle configurations, (e.g., NAR configurations 129 and 134) 
the delta portion is preceded by a cylindrical fuselage. In this region, the 
swept cylinder theory of Beckwith, and Beckwith and Gallagher should be used. 

On the wing and at low angles-of-attack, two dimensional methods are probably 
sufficient but at higher angles, cross flow effects should be included following 
the procedure of Beckwith and Cohen. 

To be conservative, it is recommended that "spot" calculations be made 
to determine upper and/or lower bounds in terms of 2D/3D and laminar/turbulent 
predictions. This of course would not be necessary if exact solutions were 
available. In addition the degree of confidence in the overall predictions 
is increased as exact methods for local predictions become available. The 
recommended procedures, which are summarized in Table 1, are not expected to be 
applicable to all possible vehicle configurations but are considered adequate 
for the current delta configurations. In using the recommended procedures 
or any other approximation scheme, it should be noted that, although a number 
of possible "improvements" are available, the increased complication is often 
not warranted in view of the fact that these improvements are building onto 
solutions that do not account for all of the important phenomena. The test of 
the accuracy of approximation methods is a comparison with exact solutions or 
flight data. Because of the budgetary and political constraints imposed on 
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shuttle, it is not practical to wait for flight data, hence the development 
of exact or near-exact computer codes (at least for specific regions of the 
shuttle surface) is the key to a successful shuttle design. 
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SECTION 5 

CANDIDATE ADVANCED BOUNDARY LAYER TECHNIQUES 

During the past few years with the advent of large scale computer machinery 
a number of computer codes have been developed which yield numerical solutions 
to the boundary layer equations. These codes are relatively expensive to operate 
and are often tempermental . For these reasons it is not envisioned that they 
should take the place of experiments or simple engineering relations for predict- 
ing beat-transfer rates when such methods are available. However, they are very 
useful in auxiliary studies where they can be used to validate or calibrate 
engineering relations, to extend engineering relations to include additional 
effects, to extrapolate wind-tunnel tests to flight conditions, and to develop 
correlations for use with simple engineering relations. A classic example 

of this last approach is the correlation of stagnation point boundary layer 

2 

solutions by Fay and Riddell. 

While boundary layer computational technology has come a long way, 

no code is available (nor would it be practical to develop a code) which treats 

the complete shuttle boundary layer precisely. This would require consideration 

of a three-dimensional boundary layer with separated flow and viscid-inviscid 

coupling over the leeward side of the vehicle. The only three-dimensional 

code available today which is oriented to flight vehicle gecmetries and which 

71 

treats the equations precisely is that of Der. However, this code uses an 
explicit approach and therefore is very expensive to operate. Also, it does 
not consider separated flow, viscid-inviscid coupling, turbulence, or chemical 
reactions . 

A rather extensive review of boundary layer computatic technology 
was presented recently in Reference 78. This report discussed in detail the 
various numerical procedures presently in use (e.g., shoot and hunt, 
quasilinearization, streamwise integration, implicit and explicit finite differ- 
ence procedures, and successive approximation and Newton-Raphson iteration 
methods) . It also discussed the current status of all then known codes for 
treating the chemical state, coupling to inviscid flow, coupling to surface 
phenomena, molecular transport, turbulent transport, and three-dimensional flows. 
While a few new codes have appeared since this report was written in late 1968 
and some of the codes discussed therein have been developed somewhat further, 
the basic conclusions of the report are not changed. Therefore, the reader 
is referred to Reference 78 for a comprehensive comparison of candidate codes. 
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Based on the results of this survey, it would appear to be practical 
to develop a code for shuttle application with the following features: 

1. Consider-; the full nonsimilar axisyrnmetric or planar boundary layer 
equations (i.e., it contains no similarity approximations). 

2. Solves the boundary layer equations "exactly" in a numerical sense 
(i.e., without the use of linearization, assumed profile shapes, etc). 

3. Considers laminar, transitional, and turbulent flows including rather 
sophisticated models for transition length and turbulent eddy viscosity 

4. Performs calculations along inviscid streamlines with approximate 
treatment of three-dimensional effects through use of small cross- 
flow theory (the axisyrnmetric analogy) . 

5. Considers regions of attached flow only. 

6. Considers nonisentropic boundary-layer edge expansions (i.e., entropy 
layer effects) . 

7. Considers detailed nonequilibrium chemistry including surface-catalyzed 
reactions or equilibrium chemistry when the situation dictates. 

8. Permits consideration of surface ablation materials through use of 
a relatively general chemistry model and general ablating wall 
boundary conditions. 

9. Performs each streamline solution around the body in a matter of a 

few minutes on a relatively high speed computer such as the Uni vac 1108. 

There are possibly several dozen separate codes in the country today 
which perform Items 1 ant? 2 above, that is, which yield "exact" numerical 
solutions to various sets of boundary layer equations. However, most of these 
are limited to incompressible or compressible single-component boundary layers 
or have some other major shortcoming that rules them out of consideration for 
application to Shuttle. Experience in developing the BLIMP code at Aerotherm 
would suggest that the development of an operational compressible single- 
component code is a trivial fraction of the effort needed to develop a code 
of the type needed for Shuttle. 

Limiting attention, then, to reacting, nonsimilar boundary layer codes, 
there are about six codes (or classes of codes) which would appear to have 
some potential for application to shuttle. These are 

1. Implicit method of Blottner as extended and used by several investigators 

2. Shoot and hunt method of Smith 

3. Implicit method of Cebeci 


36 



4. Multiple strip method of Pallone 

5. Implicit metnod of Spalding 

6. Newton-Raphson method of Kendall (BLIMP program) 

Probably the earliest nonsimilar nonequilibrium boundary layer code 
was developed by Blottner while at General Electric. In Reference 10 he ex- 
tended the implicit finite difference procedure developed while a Flugge-Lotz 
79 

student to binary nonequilibrium (atoms-molecules) . The basic computational 
approach developed under this thesis was used later by several subsequent stu- 
dents including Fannelop and Davis. The Blottner binary nonequilibrium 

code was obtained by Boeing and extended by Tong to muxticomponent air nonequil- 

82 11 
ibrium. Blottner also extended his code to a nonequilibrium tiir model. 

This code has been used and extended at ARO by Adams, Davis®^ and Lewis®"* 

in studies of chemical nonequilibrium, mass transfer, viscous interaction and 

turbulent flow (equilibrium only) . The General Electric version has also been 

86 87 

extended to include ablation products ' and to apply to the thin shock 
layer equations. 88 

Another arly nonsimilar boundary layer code which has gotten extensive 

8 9 

use was developed by Smith and Clutter of McDonnell Douglas. This code 

78 

uses a shoot and hunt method. It was extended to nonequilibrium air in 

Reference 90 and to include binary (foreign gas) injection (but for equilibrium) 

in Reference 91. This latter version has been used, for example, by Mayne 
92 

and coworkers. Further extension of these codes seems to have stopped with 
the advent of an improved implicit finite difference method by Smith and 
Cebeci. ' This code has been used rather extensively by Cebeci (e.g.. 

Ref. 95) in turbulent model studies including air-to-air injection but while 
still retaining a nonreacting compressible boundary layer framework. Other 
compressible turbulent nonreacting boundary layer codes have been developed and 

96 

used m turbulent model studies, the most noteworthy being Bushnell and Becxwith 

97 

and Herring and Mellor. 

A third early approach for solving the nonequilibrium, nonsimi’ar boundary 

98 

layer equations was developed by Pallone and covorkers. They employed a 

78 

multiple strip integral method. This code dees not appear to have been used 
too extensively, but some calculations are reported in Reference 99 considering 
teflon ablation products and further solutions are promised (but not presented) 
in Reference 100. 

In 1967 Spalding and Patankar published a book 101 containing listings 
and instructions for use of a skeletal computer program based on an implicit 
finite difference procedure for solving chemically reacting laminar or turbulent 
flows, the idea being that the user would supply the necessary subroutines 
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(chemistry, turbulent model, etc.) to solve his specific problems. This 
method has been used extensively by Professor Kays of Stanford University and 
his students in studies of the turbulent (principally incompressible) boundary 
layer. 10 ^ It was also used by Mayne and Adams 10 "* in a study of streamline 
swallowing in laminar boundary layers. No other American references could 
be found. Certainly a lot of effort would be required to develop this method 
to the point where nonequilibrium chemistry is included. 

78 

In 1966 a novel implicit procedure employing Newton-Raphson iteration 

104 

was developed by Kendall and Bartlett of Aerotherm. This method ha e been 
used extensively to study equilibrium chemically-reacting laminar ana turbulent 
boundary layers with and without surface ablation, including entropy layer 
effects, and including the axisymmetric analogy to three-dimensional flow. 

A fairly recent version of the program which is termed BLIMP s described in 
Reference 105. While the code is currently operational for general equilibrium 
flows only, subroutines governing homogeneous nonequilibrium are currently 
built into the code which, while not fully operational within BLIMP, are operational 
as a separate nonequilibrium streamtube code. 

Additional code developments worthy of mention include those of Galowin 

and Gould of GASL, 106,107 Mondrzyk of Boeing, 108 Moore and Lee of TRW, 109,110 

Chenoweth of Sandia, 111 and Marvin and Sheaf fer of NASA Ames. 86 The first two 

codes treat nonequilibrium but employ time consuming explicit procedures. They 

appear to have been used very sparingly. Moore and Lee present solutions for 

discontinuous inert injection into a nonequilibrium laminar boundary layer 

but no solutions seem to have been reported since their initial efforts. 

Chenoweth presented plans for a nonequilibrium turbulent boundary layer code 

but it is not clear that the code has ever been developed. The code of Marvin 

and Shcvffer, while limited to a binary equilibrium mixture, is mentioned 

112 27 

since it is being used to evaluate shuttle heating data. 

While the codes mentioned above consider for the rnosi: part reacting 

nonsimilar flows, they are not all ideally suited for application to shuttle. 

108 

Some apply only to sharp bodies and some require special starting procedures 

98 

(eg., marching from a known solution). Many are currently limited to laminar 

flow (e.g. , Refs 86, 99 and 110) and little experience has been gained with 

li 90 98 

some (e.g., Refs 101, 106, 110 and 111). Many consider only air chemistry ' ’ 
While this would be sufficient for nonablating heat shields, it would be 
severely limiting in the event that replacable ablation panels come into vogue. 
Th°re are several other considerations such as calculation^ speed, ease of 
setting up and running problems, and generality — in other words, usability. 
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Taking all of these factors into consideration it appears that there are two 
major candidates for a shuttle boundary layer code — the BLIMP code and a re- 
cent version of the Blottner code. 

Extension of the BLIMP code would require the implementation of the homo- 
geneous kinetics subroutine and other technical but straightforward modifications 
required to achieve nonequilibrium solutions with nominally the same accuracy, 
staLility, and computational speed with wh5ch equilibrium solutions are currently 
generated. Extension ot a Blottner code would require (depending upon the spe- 
cific starting point) the implementation of a transitional and turbulent heating 
model, addition ->f the axisymmetric analogy, and the addition of a general abla- 
tion capability. Judging from the number of codes which were developed in the 

11 90 98 

early 1960 S s to include nonequilibrium ' ' and the fact that turbulence and 

abla.lon phenomena have been included only recently, it would appear that the ex- 
tension of BLIMP to nonequilibrium would be substantially easier. Secondly, the 
resulting code would be apt to be considerably more general and flexible since 
this is the major advantage of the BLIMP code over other codes today. Finally, 
based on the experience of one of the authors (HT) who has used both codes, the 
BLIMP code would be expected to be easier to use and less expensive to operate. 
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SECTION 6 


BLIMP PREDICTION CAPABILITY 

BLIMP is a nonsimilar, chemical equilibrium computer code which uses a 
spline fit technique for approximating the distribution of flow variables within 
the boundary layer. The present capability includes an arbitrary edge pressure 
distribution, surface kinetics, laminar and turbulent flows, two-dimensional 
entropy layers and cross flow (axisymmetric analogy) . All of the above features 
have been validated either on wind tunnel models or Apollo flight data. The 
spline fit technique is very efficient in terras of computational times when 
compared with finite difference methods. 

In its present state, BLIMP can be used to predict pitchplane heating 
rates to the windward and leeward (for no separation) surfaces during the 
equilibrium chemistry portion of the trajectory. To make BLIMP applicable to 
the rest of the trajectory, nonequilibrium chemistry and three-dimensional entro- 
py layers should be included. Even without these modifications it is believed 
that nonequilibrium heating rates to catalytic surfaces are not significantly 
different from equilibrium heating rates, so that with a bit of empiricism, 

BLIMP c^uld be used for the full trajectory provided the shuttle surfaces are 
catalytic . 

For noncatalytic surfaces, BLIMP in its present state is inadequate 
and nonequilibrium chemistry of some sort is essential. Since the primary 
driving function for heat transfer is enthalpy, an elaborate multi-component 
representation is probably not required. It has been shown that binary models 
are sufficient in most cases, but for shuttle, since oxidation is important, 
a minimum model would require at least three species, namely atomic oxygen, 
atomic nitrogen, and molecular "a'r". With nonequilibrium chemistry, BLIMP 
would then be able to handle surfaces with discontinuous surface catalycity and 
oxidation rates. 

BLIMP, like other boundary layer codes, requires a means of either 
determining or specifying the edge conditions. A scheme such as that used by 
Adams can be used since BLIMP's present capability allows for streamline 
absorption into the boundary layer. Alternatively, whe methods described in 
Section 4.4 can be used. 
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Since BLIMP account? for most of the phenomena associated with entry 
vehicle heating, a favorable comparison with wind tunnel test data generates 
a high degree of confidence in flight predictions. The accuracy of these 
flight predictions have been proven for Apollc class vehicles and this experience 
suggests that a code such as BLIMP should be used in a sensitivity study and 
the results compared with wind tunnel data to determine which basic assumptions 
and phenomena will be of significant importance in flight. It should be stressed, 
though, that wind tunnel conditions will generally not include chemical dissoci- 
ation-recombination effects so that the ability of a code to adequately account 
for these effects must be determined from extensive experience with the code. 

To assess the effects of various three-dimensional and entropy layer 
assumptions, a matrix of heat transfer solutions were obtained for the windward 
generator of the NAR Model 134B and compared with the wind tunnel data of 
Reference 27 and 113. Additional supporting data were obtained by direct 
communications with the authors of Reference 27. The body geometry of the NAR 
Model 134B was obtained from Reference 114 and shock shape was measured from 
shadowgraphs of the flow about NAR Model 129^'* (which is slightly larger than 
Model 134B) . These shadowgraphs were obtained from J. Cleary of NASA Ames, 

Since the wind tunnel test conditions were at relatively low temp, ‘■ures, 
the homogeneous version of BLIMP which consumes substantially less compiler 
time than the chemically reacting version was used. These short run times 
make homogeneous BLIMP an economically practical tool for extensive studies 
of wind tunnel test conditions. Several check runs wer made with the chemical 
code to verify the accuracy of the homogeneous results and in all cases, 
predicted heat fluxes were found to agree within ± 2%. 


6.1 PRESSURE DATA FOR BLIMP INPUT 

Pressure ratio data (P/Pt 2 ) were taken directly from the pressure ratio 
plots in Reference 113. The total pressure behind the shock was obtained assuming 
an ideal gas with y = 1.4 and M = 7.4. It is noted that there is undoubtedly 
some error introduced by interpolating from such a plot; however, the uncertainty 
introduced into the heat flux predictions will be only about 50 percent of the 
uncertainty in pressure and thus this source of error is considered small (no 
more than 3 percent at the lowest angle of attack) . For the purposes of the 
BLIMP code input, additional pressure data were required, particularly in the 
nose regions. These data were generated by assuming a Newtonian pressure 
distribution from the stagnation point and calculating the pressure ratio from 
the measured slope of the windward aenerator (see Section 6.1.2). The pressure 
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distributions used in the predictions are given in Appendix 2 and consisted of 
these computed values faired into the forwardmost reported measured values. 
Pressure gradient data (3(P/P T )/3s) for the spreading factor calculation were 
obtained by measuring slopes from the pressure ratio vs. s plots. 

Heating data were obtained from the normalized plots in Reference 113. 
in a similar manner. Actual heating values were computed using the calculated 
stagnation point heat flux values listed in Table 2. These values were computed 
by Marvin based on a 0.006 foot radius sphere and were used to normalize his 
data. These values are not measured stagnation fluxes. 

6 . 2 MODEL GEOMETRY 

The side profile and several frontal cross sections of the NAR Model 134B 
are shown in Figure 8 to provide a general indication of the model configuration. 
Pertinent geometric variables are 

s , distance along the windward generator at each angle of attack. 

Here it was assumed the stagnation point was located at the inter- 
section of the windward generator and a radial line through the 
center of the nose sphere at the particular angle of attack. The 
body centerline was taken from the drawing co be a line through 
the center of the nose sphere and parallel to the straight sections 
of the upper and lower fuselage. 

R c , radius of curvature of the windward generator 

R t , radius of curvature at the windward generator in the plane transverse 
to the generator. This plane is locally perpendicular to the 
generator . 

r Q t the perpendicular distance from the windward generator to the 

angle of attack axis, i.e., the local radius of an axisymmetric 
body generated by rotation of the windward generator about the 
angle of attack axis. 

6.3 SHOCK WAVE GEOMETRY FOR ENTROPY LAfER PREDICTIONS 

Shock wave geometry was obtained from full scale shadowgraphs of Model 129 
which were obtained at the same freestream conditions as used in heat transfer 
and pressure tests. Angle-of-attack conditions included a = 15°, 30°, 45° 
and 60*. Data for the o = 53.5° test condition were obtained from interpolation 
between the a = 45° and 60° cases. 
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The shock angle (given in Appendix 2) relative to the angle of attack 
axis was measured as a function of the radial distance from the angle of attack 
axis to the shock and then converted into total pressure ratio across an oblique 
shock by the relation 



(y + l)M^sin 2 0 

Y-l 

y + 1 

(y - l)M^sin z 0 + 2 


2yM^sin z 0 - (y - 1) 


(6-1) 


This ratio together with the shock radius was used for the entropy layer input 
for the homogeneous code. 


6.4 APPROXIMATION OF THREE-DIMENSIONAL FLOW EFFECTS 

Three dimensional effects can be approximated in the BLIMP code by an 

• * 

axisymmetric analogy or in terms of a streamline spreading factor h 2 . In 
the former approximation, the heat transfer rates are calculated by assuming 
that the body is axisymmetric with a longitudinal profile specified as the 
windward generator. The continuity equation is then given by 


7S <h 2 P u > + ^ <h 2 Pw> 


0 


For a spherical surface such as on the face of an Apollo vehicle, fc 2 is given 11 ** 
by the equation 



1 3P ^2 . 
pu 2 3s ds 




( 6 - 2 ) 


where 


y B is the angle between the local surface normal and the free 
stream velocity vector 

P T is the total pressure 

P- c is the local surface radius of curvature 


For this case, there is no ambiguity in the variable R since the surface is 
spherical, but for a general chree-dimensional body such as a shuttle vehicle 


h 2 is also the metric coefficient of the streamwise coordinate 
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the above formulation does not account for the difference in surface curvature 
as measured in the transverse and longitudinal directions A new derivation, 
presented in Appendix 3, shows that for non-spherical surfaces the spreading 
factor should really be defined by the equation 



1 

pu 2 


ap _ 2. 

3s ds 



2P t cos 2 0 b 



0 


(6-3) 


Equation (6-3) was integrated using a Runge-Kutta method 
second order equation. The dynamic pressure term was evaluated 
gas relation 



for a single 
from the perfect 


(6-4) 


resulting in the equation being a function of the two surface radii of curvature, 
the pressure ratio, and the pressure ratio gradient which was evaluated from a 
plot of the measured pressure supplemented by the calculated Newtonian pressure 
gradient distribution near the nose. Calculated values of hj are given in 
Appendix 2. 


6.5 BLIMP PREDICTION MATRIX 

Experimental data were available for laminar flow at a = 15°, 30° and 
53.5° and turbulent flow at a = 30°; the matrix of BLIMP solutions, shown in 
Table 3 , is limited tc these angles-of-attack and the test conditions shown in 
Table 2 . 

6.6 RESULTS AND DISCUSSION OF BLIMP STUDIES 

Figures 9 through 12 show the influence of various assumptions on the 
predicted heating rates. The small crossflow results shown are all based upon 
P/P T replacing cos 2 0g in the spreading factor equation (Equation (6-3)) and 
R^ = » in the aft or wing region of the fuselage. All predictions for a given 
case used the same stations, pressure ratios, and wall temperatures input distri 
butions. Three-point differencing was used throughout to obtain streamwise de- 
rivatives except for the turbulent case which used two-point differencing.* 


A run using three-point differencing was also made with negligible difference 
in the results in the laminar region. Two-point differencing was used to provide 
more flexibility through the transition region than afforded by the limit of 
4, two-point difference stations in the three-point difference option. 
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Laminar Flow Results 


The blunt planar body predictions made for the a = 15° and 30° (Figures 
9 and 10) laminar cases under predict heating over the entire length of the wind- 
ward generator. The differences increase with angle of attack to roughly 
50% for a = 30°. These trends are consistent with the over-prediction of 
boundary layer growth resulting from the planar approximation as compared to 
an axisymmetric treatment. Although the planar blunt body assumption may be 
useful in establishing a lower limit on the predicted heat transfer rates, its 
use is not justified since there is no significant reduction in time to either 
set-up the input data or run the code. However, since no wing heat transfer 
predictions were attempted, no conclusions can be reached on the validity of 
this assumption in these regions. 

The predictions based on an isentropic expansion around a blunt axi- 
symmetric body compare favorably with the data for the a = 15° and 30° 
cases (Figures 9 and 10), particularly over the aft regions of the fuselage. 

In the nose and forward regions of the fuselage, the predictions are somewhat 
low; an apparent effect of the lower degree of streamline spreading caused 
by the assumed axisymmetric body as compared to that which should exist over 
the smaller transverse curvature of the actual body. This lower spreading of 
the flow causes a more rapid boundary layer growth and consequently reduced 
heating. This effect is confirmed and amplified in the a = 53.5° case (Figure 
11) where the predictions are about 50% low in the forward region. In addition, 
it appears that the history of this thicker boundary layer also causes an 
underprediction of the heating along the aft regions. 

The results of the third assumption, that of the small crossflow theory, 
show excellent agreement with data in the forward regions for all three angles 
of attack. There is a tendency to overpredict wing region heating at a = 15°. 

At a = 30°, only the first half of the wing region shows any disagreement. 
Finally, at a = 53.5° the small crossflow prediction is low along the latter 
half of the body by an average of 25*. It is noted that Marvin, et. al., 
(Reference 113) suspect that the final data point at a = 53.5° (flagged in 
Figure 11) is in a region of transition to turbulent flow. BLIMP-predicted 
values of Reg at this location were of the order of 200 to 250. Based on the 
correlations given by Masek 6 ^ for the dependence of the transition parameter 
Reg/[M e » (Re/L) ] 0 • ^ with angle of attack, transition might be expected at Re Q 
of approximately 150 to 175. If the flow were laminar, the heat flux at this 
last body location would be less than that of the previous upstream location. 

For a blunt axisymmetric body assumption, accounting for an entropy layer 
(as opposed to an isentropic expansion) results in an increase in the predicted 
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heating rates of about 10%, 15% -and 25% respectively for the « = 15°, 30° 
and 53.5° cases. Inclusion of an entropy layer appears to be only a small 
improvement over the isentropic expansion case, but even these small effects 
may have significant influences on a shuttle vehicle design. A measure of 

the rate at which flow from the high entropy stagnation region shock is 

swallowed by the boundary layer is obtained by comparing the shock wave angle 
through which the boundary layer edge streamline originated. For a = 30%, this 
shock wave angle is plotted in Figure 13 as a function of the surface coordinate s 
and it can be observed that the swallowing process is essentially complete before 

s *v» 0.1 ft. The continuing but slow decrease of the shock angle, and thus the 

entropy, is primarily due to the continual decrease of the shock angle along 
the entire body length (see Appendix 2) . Figures 14 and 15 show typical boundary 
layer velocity profiles (solid lines) for a = 15° and 30° at x/L of 0.7 to 0.8. 

The dashed curve represents the shock angle of the streamline as a function 
of y, the normal distance from the wall. These plots also indicate that the 
entropy layer has been swallowed well forward on the body. 

Predictions combining the entropy layer with the small-cross flow theory 
were not made as this system. has not been incorporated into the current version 
of the BLIMP code (this is conceptually straightforward but requires some code 
modification). However, by extrapolating the available results, small-crossflow 
with entropy layer will slightly increase the isentropic small-crossflow results 
for a - 15® and 30° and should substantially improve the prediction over the 
latter half of the body for the a = 53.5' case. 

Turbulent Flow Results 

The turbulent flow results for a = 30° are presented in Figure 11. 
Transition to turbulent flow in the wind tunnel experiments of Reference 113 was 
reported to occur at x/L = 0.54. Transition to the fully turbulent model in 
the BLIMP code was made to occur at the ;ame location, thus the apparent 
instantaneous increase in turbulent heating levels. The interesting result 
here is the effect of the entropy layer on turbulent heating. Whereas in the 
laminar region the increase, above the isentropic axisymmetric value, is limited 
to about 10%, in the turbulent region the increase is 40% resulting in good agree- 
ment with the data downstream of the actual transition region. This difference 
in relative effects is most probably due to the greater dependence of transport 
properties near the wall on edge conditions in turbulent flow as compared to 
laminar flow. At the transition station edge gradients still persist primarily 
due to shock curvature along the body .is shown by Figure 16 which includes the 
laminar velocity profile and the streamline shock tagle at the final laminar 
station. The reported transition momentum thickness Reynolds number. Re , in 
Reference 113 was approximately 450. Predicted values for the fixed transition 
location were 470 for the axisymmetric rur. and 330 for the small crossflow theory. 
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OTHER RESULTS 


Equations (6-2) and (6-3) were derived with a Newtonian pressure assump- 
tion. For shuttle configurations, this is not quite correct so that the validity 
of these equations might be questioned. In Equation (6-3) , the Newtonian 
assumption appears solely in the second term of the hj coefficient where P T cos 2 0 B 
represents the Newtonian pressure at's on the plane of symmetry. At a given 
0g the actual pressure will be slightly higher; in particular, as 0 fi approaches 
90°, P t cos 2 0 b goes to zero whereas the real pressure does not. Consequently, 
replacing the term P T cos 2 0 B with the actual surface pressure would increase 
the magnitude of the corresponding term in Equation (6-3) . 

To investigate this point further, two solutions were obtained for the 
a = 50° spreading factor case, one with an distribution based on P T cos 2 0 B and 
the other on the measured pressure on the symmetry axis. The two results were 
identical from the stagnation point to an x/L of 0.1. From this point to the 
beginning of the fla at x/L = 0.2 the heat fluxes based on the measured pressure 
increased to 5 percent above the P T cos 2 0 B value. Over the remaining (wing) 
portion of the body the measured values were from 0 to 2 percent higher than 
P T cos 2 0 B . Thus, it is concluded that for the shuttle configuration, this effect 
is of a minor order. For all angles of attack, the reported results are 
based on use of the local pressure which should overpredict local heating by 
no more than 2 percent. 

The effect of transverse curvature in the aft region was evaluated by 
assuming that R^ was (a) infinite and (b) equal to that of a cone locally 
tangent to the surface and coaxial with the angle-of-attack axis. For a = 30°, 
the tangent cone approximation resulted in an increase in heat fluxes from 5-10%. 

Additional comparisons which would be beneficial are shown in Figures 17, 

18 and 19. In these figures the BLIMP axisymmetric blunt body predictions are 
compared with the tangent cone approximations recommended in Section 4.4. The 
data of Reference 113 are also included. At all three angles-of-attack (15°, 

30°, 53.5°) the tangent cone approximations are slightly higher than the BLIMP 
prediction but both are in generally good mutual agreement. However, at 
a = 53.5° both predictions are noticeably lower than the data and would suggest 
that other factors such as three dimensional entropy swallowing or variable 
transverse curvature effects are important, especially at high angles-of-attack. 
These effects are presently approximated in the BLIMP code. Although these 
approximations do not necessarily represent optimum approaches, improved 
methods will probably rely on the development of near-exact inviscid solutions 
such as the work of References 21 and 22. Consequently an important consideration 
in the development of a boundary layer code is the inclusion of an option for 
integral coupling with an inviscid code. 
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SECTION 7 


RECOMMENDATIONS FOR FURTHER DEVELOPMENT OF 
HEATING PREDICTION CAPABILITY FOR SHUTTLE VEHICLES 

Although it is possible *-.o predict the approximate magnitude of the heat 
transfer to be expected on a delta shuttle vehicle, the reusability concept 
dictates a higher level of accuracy. The following are some recommendations 
on methods and procedures to improve accuracies of present methods. 

1. The BLIMP code should be modified to account for nonequilibrium 
homogeneous chemistry. 

2. Nonequr librium BLIMP should be used to evaluate the influence of 
locally noncatalytic surfaces on downstream heat transfer. 

3. The BLIMP code should be used to evaluate the influence of various 
approximations of streamline spreading and entropy swallowing 
effects for typical shuttle flight conditions. 

4. A critical examination and evaluation of transition-to-turbulence 
information should be performed to define a viable transition criteria. 

The present BLIMP code has options for integral coupling with inviscid flow 
field calculations and surface reaction chemistry. Any modifications and 
improvements to this code should be developed so as to retain these options. 
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APPENDIX 1 

LOCAL SURFACE INCIDENCE RELATIVE TO 
FREE STREAM DIRECTION 

Local Surface Incidence of Wing Surface 
Defining pertinent angle as 

a = vehicle reference angle of attack 
a = angle oetween chord line and reference line 
a = local surface inclination with respect to chord line 
$ = local dihedral angle 



then at a = 0, the unit vector which is tangent to the surface and ne s in the 
chord section plane is 

A ^ = A _ A 

= cos (a + a)i - sin a j 


and the unit vector which is tangent to the surface but which lies in the plane 
perpendicular to the reference line is 


I 

I 


Ug = sin <j> j + cos <f> k 


tfl* 
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The cross product of these two vectors is normal to the local tangent plane and 
in normalized form is 

f * _ sin (a + a)cos 4> ? ccs (a + a) cos <j> ? ^ cos (a + a) sin <J> 

U 3 = N 1 N 3 N 

where 

N 2 - sin 2 (a +a)cos 2 <f> + cos 2 (a + a)cos 2 <f> + cos 2 (a + a)sin 2 if> 

=s cos 2 4> + cos 2 (a + a)sin 2 <f> 

\ 

Rotating through an angle of attack a 



A 

The length of the projection of onto the x-y plane is 


and 


COS (j> 
N 


= a + a 
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The unit normal vector which defines the local tangent plane at angle of attack 
a is then 


* _ sin (a + a + a) cos <p T cos (a + a + a) cos <t> * , cos (a + u + a) sin it> 

u 4 N 1 N 3 + ~ 

N 

A A 

The dot product of with i is the sine of the effective xncidence. Defining 
Y as the effective incidence 


sin y = 


_ cos <{> sin(a + a a) 


N 


where the sign was 
For small a 


changed to yield a positive angle. 
+ a this reduces to 


sin y “ cos ^ sin (a + a + a) 


Local Surface Incidence of Outboard Fin Surface 
Defining pertinent angles as 

a = angle between fin chord line and vehicle symmetry plane 
= angle of tilt of fin 

a' = local profile angle relative to chord line 

then at a = 0, the unit vector at an angle from the vertical and in a plane 

perpendicular to the fin chord line is 

A A A A . 

= - sin 4>»sin a i + cos <j>'j + sin <J>'cos a k 

The unit vector which is tangent to the fin surface and lies in a horizontal 
plane is 

A A. A 

U 2 = cos (a + a')i + sin (a + d')k 


Vf > 



The cross product of and after normalization yields the unit normal 
vector which defines the tangent plane. Thus 


U 3 * ' 


cos ♦ 




sin a sin (a + a 1 ) 


a')| j 


+ cos ' cos (a + a* ) 


+ sin ♦ ' cos a cos (a + a') > j + — 7 k 

N 1 


where 


•{ 


N? = cos ij)' + sin 2 ^' {sin a sin (a + o') + cos a cos (a + a') 


Using the same procedure as before, is rotated through an angle a. The 
unit normal vector which defines the local surface tangent at angle of attack 
a is then 


°4 = - 


L rin (a + O i - L cos (a + 5)j + P° s -t l- g . ° s iS- y 


where 


L 2 = cos 2 < fr' s i n2 (a + a ' ) + J- 


N i 


N 


sin ♦’sin a sin (a + a 1 ) 


+ sin ♦ 'cos a cos (a + 


a') | 


tan K = 




cos ip 'sin (a + a') 


sin ♦'sin a sin (a + a’) + sin ♦'cos a cos (a + 


a’) 


Again, the dot product of and i is the sine of the effective angle of incidence 
y'. With the sign changed to yield a positive angle, 

sin y* = L sin (a + £) 


and L and £ can be expressed as 


L = 


cos^ 'sin* - (a + a*) + sin 2 ^ 'cos*a ' 


cos 2 4 { + sin 2 ^ 'cos 2 (a* ) 


^ A * n/ie ^ t r* * 


1 

*> 


tan 


-l 


sin (a + g'l 
cos a* tan ♦' 
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I 

L 


! 0 
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i 
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5 


t li 


i 

f- 


For small angle a' this set of equations reduce to 


sin y' = u 


s <j> ' y sin 2 a + tan 2 <J> 1 sin(£ + a) 


sin K - 


sm a 


^sin 2 a + tan 2 4>' 


True Sweep Angle of Wing Leading Edge 
Define 

t = semi-apex angle of wing 
$ = dihedral angle 

Then consider the unit vector which is the axis of a hypothetical cylindrical 
leading edge. This unit vector is 

r, _ cos C 7 * sin C tan 4> ? ^ sin C C 
U 1 " N i ~N^ 3 N x K 


where 


= 1 + sin 2 ^ tan 2 $ 

Rotating this vector through an angle of attack a, the unit vector becomes 


where 


U 0 = L costa - A,)i - L sin (a - A~)j + 


L 2 - cos 2 C + sin 2 C tan 2 <{) 

N? 


sin c 
N, 


1 I 


1 B 

xx 


I 


tan A 9 = tan C tan $ 

r> A 

the dot product of U 2 and i is the cosine of the effective angle of attack A 

cos A = L costa - Aj) 

For small angles $ this reduces to 

cos A - cos C cos a 


Thus 
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True Sweep Angle of Outboard Fins 
Define 


fi = nominal sweep angle of fins measured or. vehicle symmetry plane 

also <J>' and a are as previously defined. Then the unit vector which represents 
the axis of a hypothetical cylindrical leading edge is 


£ _ sin 6 7 . cos 5 7 . (sin 6 tan a + cos 6 tan ' ) C 

U 1 N x N 1 3 **1 

where 


N 2 = 

Rotating through an angle 

U 0 = sin (a + 6)i 
2 N 1 


1 + (sin 6 tan a 

of attack a, the 

1 * 
+ rr- COS (a + 6) j 

“l 


+ cos 6 tan 4 > ' ) 2 
unit vector becomes 


sin 6 tan a + cos 6 tan <j>' 



and the effective angle of attack of the fin leading edge is 


cos A* 


sin(a + 6) 

^1 + (sin 6 tan a + cos 6 tan <j> ' ) 2 
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APPENDIX 3 


SPREADING FACTOR FOR NONEQUAL BODY CURVATURE 
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SPREADING FACTOR TOR NONEQUAL BODY CURVATURE 
The variable R in Equation (6-2) is not ambiguous since it applies only 

V 

to a sphere, the application being the Apollo reentry heat shield. For shuttle 
application, Equation (6-2) is not valid since it does not account for the un- 
equal surface curvatures which exist forward of the wing region. The deriva- 
tion of Reference 116 assumes Newtonian flow along a streamline and transverse 
to it in arriving at an expression for the crossflow pressure gradient in the 
meridianal direction. 


3P 

3m 


COS ! «P 

K 


(A-l) 


where z is the distance from the plane of symmetry. 

For a non-spherical surface, it is necessary to return to the basic 
assumption of Newtonian flow, namely 

P(s,z) _ cos 2 ip (s ,z) . 

P(s,0) cos 2 iMs,0) 

where ip is the angle between the local surface normal and the free 

stream velocity vector. 

Referring to the sketch, the circle with radius R^ represents 


WiMPWAEP 
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the tangent circle in the plane perpendicular to the windward streamline 
(generator), z = 0. A rotation in this plane by an angle g, which is equiva- 
lent to some displacement z from the plane of symmetry, results in a change 
in the angle ip. This angle may be found by viewing the tangent circle both 
in planform and on edge. The point Q represents the tip of the radius 



vector at some angle g. The angle this vector makes with the plane perpendicu- 
lar to the velocity vector V is simply 


a = sm 


- 1 d 




(A- 3) 


and thus 


ip ( 3 ) =90 - a = 90 - sin 




(A- 4) 


but from the above sketches we have the relation 


d _ 
iyzosiMO) 


(A- 5) 


and 


t = R T (1 - cosg) 


(A- 6} 


thus 


i MS) = 90° - sin 


in 1 rJj, j^rC 


1 - (1 - cosg) ) > cosp (0) 


) | cosp (0) 


(A- 7) 


or 


4»(g) = 90° - sin -1 £cosg cos<|>(0) 


0 


(A- 7a) 
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Returning to the Newtonian pressure expression, Equation (A-2) 
possible to find the pressure gradient in the direction normal tc the 
streamline. 


■ 


or 


where 


3_ 

P(s,z) 

= 3_ 

COS 2 (s , z) 

3m 

P (s /0) 

3m 

cos 2 /]/ (s ,0) 



cos 


P(s,0) 

(s ,0) 


a_ 

3 m 


£cos 2 4/ (x,z)J 



sin 


- 1 


^cos3 cosij/ (0) 



but 


cos (90° - y) = siny 


and . sin (sin” l x) = x 

therefore 


(cos 2 ///) = £cos 2 3cos 2 /// (0)J = COS 2 l|) 


(0) ^ (cos 2 3) 


Now for small z and thus 3, Equation (14) becomes in the limit 


^ (cos 2 ///) = cos 2 \M0) 


[- « s] 


with 


g « and m # z 

*T 


thus 


|j (COS 2 l|0 -co.'*«0)|-2^gg 


= COS 2 4»(0) 


[■ 

(•s) 


] 


it is now 
symmetry 

(A-8) 

(A-8a) 

(A- 9) 

(A-10) 

(A-ll' 

(A-12) 

(A-13) 
j (A- 14 ) 

/ (A- 15)' 


in the limit as z goes to zero 



’'•vV 




Jf the lojal symmetry plane pressure in Equation (A-8a) is the Newtonian 
value then 


P(0,0) = P 

cos 2 ^ (s,0) 


and Equation (10a) becomes 


~ cos 2 \l> (0)P T 
R T 


where (0 ) and 0,- are equivalent quantities. 

Thus, Equation (6-2) is replaced for non-sphetical shapes by 


1 3P ^2 
...2 3s ds 


" j, _ !!i cos26 b 1 h = 

y c pu 2 r 2 J 2 
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